Quantum entanglement of flux qubits via a resonator. 
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Abstract 

We show that flux qubits can be efficiently entangled by inductive coupling to 
a tunable resonant circuit, in the scheme reminiscent of atoms' entanglement 
through the optical cavity mode. It is shown, in particular, that the single- 
photon excitation of the resonator produces the pure Bell state of qubits with 
the completely disentangled LC circuit. 
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1. Successful realization of quantum computing critically depends on our ability to 
entangle quantum states of qubits fl]]. Controllable coupling mechanism necessary to perform 
this task must also allow for variation of single qubit's parameters, which is inevitable if 
qubits are macro- or at least mesoscopic devices. This is the case for superconducting flux 
qubits [0,0. The alternative approach, fine tuning of single qubit parameters, is not always 
applicable. 

The natural way of coupling superconducting qubits is through an intermediate resonant 
LC circuit (tank) 0-0. The tank circuit can be made tunable by including in it a biased 
Josephson junction 0,0, the effective inductance of which depends on the bias current. 

In the limit 0, when the tank resonant frequency ujt ^ A, where A is a characteristic 
frequency of the qubit, the role of the tank is reduced to creating a direct qubit-qubit 
interaction, which can only be controlled by tuning individual qubits. Here we consider the 
case when ut and A are comparable, so that there are actual excitations (photons) in the 
tank. 

The presence of actual excitations in the LC contour allows more freedom, e.g. its 
consecutive coupling to the qubits, when the frequencies of the latter differ. The problem we 
approach is analogous to the one extensively studied in quantum optics PHT0|,|5| with respect 



to entanglement of atoms through the cavity modes. We concentrate here on developing 
the theory of dynamic evolution of the "two qubits plus resonator" system in the case of 
simultaneous coupling the qubits to the resonant LC circuit. The goal is to find the conditions 
which allow to maximally entangle the qubits with each other, while disentangle them from 
the quantum state of the tank. The decoherence effects are not taken into account explicitly. 
We must assume anyway that the decoherence time of qubits significantly exceeds the time 
necessary for manipulations with them (which is a realistic assumption in view of recent 
experimental successes with superconducting qubits ||TTHl3|)- The LC circuit itself should 
have the decoherence time of the same order as qubits. 

2. We consider two phase qubits, "a" and "6", having the same tunneling splitting. Each 
of these qubits is inductively coupled to the LC circuit ("tank") tuned in resonance with 
the qubits. The Hamiltonian of the tank with an inductance Lt and a capacitance Ct has 
the form 

where qt and It are a charge and a current in the tank, respectively. It should be noted, 
that for these variables we have the commutation rule: [Jt,5't]- = ifi/Lx, which follows 
from the relation between a flux in the tank, $t = LtIt, and the charge qt '■ [^t, Qt]- = i^- 
In a two-state approximation a persistent current in the qubit loop takes two values ±/a 
(for a— qubit) and can be described by the Pauli matrix ax- As a result, the operator of the 
coupling energy between the qubits and the tank will looks like: 

Hir^t = -{MJaOx + M,hTx)lT. (2) 



Here Ma = kayLoLr stands for the mutual inductance of the tank and the a— qubit having 
the inductance La- The same is true for the b-qubit which has the inductance Li, and is 
describes by the Pauli matrices t^, Ty, r^. The qubit currents la, h are constant parameters 







here. The dimensionless coefficient k obeys a relation: k'^Q ^0.1 — 1 with Q being the 
quahty factor of the LC-circuit. 

The tank represents itself a single mode cavity with excitations which we will call there- 
after "photons". Because of this, it is convenient to work with creation-annihilation opera- 
tors of photons in the tank b, where: 



V 2hujT 

Then, for the electric current in the tank we have 



'iT + ioJTqr)- (3) 



Ir-{^M. (4) 

The Hamiltonian of two qubits with the same tunneling elements = A and A^ = A 
interacting with the tank can be written as 

H = ^a, + At, - (A„a^ + XbT^){b + b+) + hurib+b + 1/2), (5) 



where the coupling constants are = Malay ^I'OOt/'^Lt, Xb = Mbh^Jii^T/'^L'r- Without the 
interaction between the qubits and the tank the evolution of operators a± = i: icTy, t± = 



Tx ±iTy,b is described as 



a±{t) = exp(±2iAt)a±(0); r±(t) = exp(±2iAt)r±(0); b{t) = e-"^^%{0). (6) 

To analyze the case of strong resonant coupling between the qubits and the tank, when 
the frequency of the tank is exactly equal to the energy splitting between the levels of the 
qubits, ujt = 2A, we use a rotating wave approximation and neglect of the fast oscillating 



terms in the Hamiltonian (5). As a result, we get the Jaynes-Cummings model [14| with 
two two-level systems resonantly coupled to the single mode cavity. The Hamiltonian of this 
model has the form: H = Hq + Hint, with 

Ho = + ^yr, + Luribn + 1/2), 

H.nt = -y + ^-b^) - y + r_6+). (7) 

As in the original Jaynes-Cummings model, two parts of the Hamiltonian H commute: 
[Hq, Hint] = 0. It means that the operators Hq and Hint have the common system of eigen- 
functions. Let's introduce aa, Pa as eigenfunction of the a,— matrix:(T^aa = aa,(JzPa = 
—Pa', and ab, Pb as eigenf unctions of the z-matrix of the second qubit, r,, with eigen- 
values ±1. Then, the action of operators a± on these functions is described as follows: 
cr+aa = 0, cr+Pa = 2aa; CT-da = 2/3a, (J-Pa = 0. The similar is true for the second qubit. The 
states of the LC circuit are represented by the eigenfunctions of the harmonic oscillator, 
\n), corresponding to the eigen-energy ujT{n + 1/2). The action of the creation-annihilation 
operators on these states is well-known: b\n) = ^/n\n — 1) , b~^ \n) = + l\n + l). We choose 
the following eigenstates of the Hamiltonian Hq (7) corresponding to the same eigenenergy 
nuT, Hq^ j = nojT^ j, i = 1, .A: 



*i = |n - l)aaab, ^4 = |n + l)f3aPb, 

*2 = \n)aaPb: *3 = HPa^b- (8) 

Here the number n is counted from 1. In the case when n — the wave function of the first 
state is assumed to be zero, = 0, and we have only three basis states. If n = 1, then in 

the first state we have no photons, but both qubits are in their upper states contributing 
(0.5 + 0.5)cjt = to the total energy (we do not consider here the tank-qubits interaction). 
In the second state, ^^2, the second qubit turns over, from the state with the energy cut/S 
to the state Pb with the energy —ur/'i. The released energy goes into the tank creating 
one photon having the energy cut- The same process occurs in going from state 1 to state 3 
when the first qubit turns over. Due to the commutativity of Hq and Hint the action of the 
coupling Hamiltonian Hint on the eigenfunctions (8) of the free Hamiltonian Hq does not 
fall outside the basis ^'i, ..^4 : 

^mt*2 = -{XaVn+ 1^4 + AfeV^^i), 

Hint'^a = -(Aav^^i + AbVn+ 1*4), 

Hint'^A = -Vn + 1(A„*2 + Ab^a) , (9) 

or, that is the same, the nonzero matrix elements of the operator Hi„f in the basis ^^i, ..^4 
will look like: 

{l\H,^t\2) = {2\Hi^t\l) = -V^K, {l\H,nt\5) = {3\Hrru\l) = -v^A,; 
(2|if,„i|4) = {A\H,nt\2) = -y^TTlA^; (3|if,„i|4) = {A\H,nt\5) = -V^TTlAfe. (10) 

The eigenfunctions {^} of the total Hamiltonian H = Hq + Hint, H'^ = E^^, can be rep- 
resented as a superposition of the basis function ^'i, ..^'4 : ^' = J^'jZtcj'^j. As a result we 
obtain the system of four linear equations for the coefficients ci, ..C4. After some transfor- 
mations these equations take the form: 

{E - uutYi = -^/n{XbC2 + AaCs), 
{E -mLiT)c4 = -^/nTl{XaC2 + XbCs), (11) 

{{E - nurf - Xln - Xlin + l)}c2 - A,A,(2n + l)c3 = 0, 

{{E - vujt? - Xln - Xl{n + l)}c^ - A,Afe(2n + l)c2 = 0. (12) 

This system has nontrivial solutions when the the corresponding determinant is equal to 
zero: 

Det ^{E- ucutT - (Aa + A6)(2n + 1){E - morf + (A^ - Xlfn{n + 1) = 0. (13) 
This condition gives us a set of four eigen-energies of the total Hamiltonian H 

E('''^ = ncoT ± ^ (Xl + XDin + 1/2) + ^(A^ + Xl^in + 1/2)^ - (A^ - X^Mn + 1), 

^ruvT±^{Xl + Xl){n + 1/2) - ^(A^ + XlY{n + 1/2)^ - (A^ - Xiyn{n + 1). (14) 

A 



3. Let's proceed to a consideration of some particular examples. In the case when the 
interaction of the tank with the second qubit is absent (A5 = 0) for the energy levels E^"^^ 
and for the corresponding eigenstates '^^'^\ = I, .A, we obtain, respectively: 



^(1.2) 



V2 



(15) 



It is evident, that the second qubit is completely disentangled with the tank as well as with the 
first qubit. Nevertheless, the first qubit and the tank are strongly entangled in this situation. 

It can be expected that the maximal entanglement could be r eached whe n both qubits 
are coupled to the tank with the same strength, Aa = = A = where Lq is a 

qubit inductance, and Iq is a current in a qubit's loop. Then the first and the second energy 
levels are different. 



= ± 2A^n + l/2, 
whereas the third and the fourth states have the same energies 

The corresponding orthogonal wave functions take the form 
= 1 J 'l—\n-l)aaab 



(16) 



(17) 



n + 1/2 

I n+1 
'n + 1/2 



I n + 1 
n + 1/2 



\n - l)aaab 



n 



n + 1/2 



n + l)/3aPb + \n)aa/3b T \n)l3aab 



\n + l)f3aPb ± \n)aaPb =F \n)(3aC(b 



[18) 



Apparently, now we have strongly entangled qubits and the tank. However, if the number 
of photons in the LC-circuit is very high, n ^ 1, so that jriil) — 1?^), the system is actually 
in a factorized state : 



(19) 



as expected, since there should be no entanglement through a classical system. 

4. It is of interest to trace the evolution of the originally non-entangled state, say, 
\l/(0) = ^0 = \0)oiaC(b- In this initial state both qubits are in the up-state and there arc no 
photons inside the tank. The time evolution of the wave function ^'(t) in the Schrodinger 
picture is described by the superposition: 



v=A 00 

*W = E E^-exp(-i£;Mi)*M, 

1/=! n=l 



(20) 



where are the eigen-energies (16), (17) and eigenstates (18) of our system. These 

states and energies depend on the number of photons in the tank n. Because of this we have 



a sum over n in Eq.(20). The coefficients Ui, in the superposition (20) are determined by the 
relation 



(21) 



and also depend on the number of photons. For our choice of the initial state the quantum 
number it should be: n = 1, and Ui = U2 = l/\/6, M3 = M4 = l/\/3. As a result the following 
four states participate in the evolution of our wave function ^(i) : 







^|2)/3aA- 


^\l){aaPb + Pa^b), 




1 , . 


■^\2)PaPb + 


^\l){aaPb + PaOib), 




1 , . 




^\l){aaPb - Pa^b), 


= 


-j=\0)aaab - 


^|2)/3a/36- 


-\l){aaPb - PaOib), 



(22) 



with the energies counted from the vacuum energy ut/^,: E^^^ — ujt -\- A-\/6, — ujt — 
Ay^, E'^^) = E^^'^ = lut, respectively. After the time interval t our system "two qubits + 
tank" will be in the state 



-^e-^'^^* sin( v^Ai) 1 1) {aaPb + PaC^b) ■ 



(23) 



It is clear from this formula that amplitudes to find the system in one or another state 
oscillate in time with the frequency 



(24) 



where Iq is the current in the qubit's loop, and k is the coupling coefficient between the 
qubit loop (with an inductance Lg) and the tank (with an inductance Lt), so that the 



mutual inductance one of the qubits and the tank, M, is: M — LqLT- The initial state, 
^'o = |0)Q;aQ;6 revives in the moments t — 27r'm/flR,m = 0,1,2,.. In the moments t — 
(7r/2 + 7rm)/Qij all states of the qubits are mixed: 



-\0)aaab- —\2)f3aPb- [-1) ^ 



(25) 



After a 1/2-period shift from the origin (at t = (tt + 2Tim) /VLr) we have an entanglement 
of two up-states of the qubits, aaO-h, two down-states, PaPbi with the vacuum state |0) and 
with the two-photon state |2) of the tank: 



1 2^2 
-|0)a„afo —\2)(3a(3b 



(26) 



It should be noted that in the states (25), (26) the qubits are entangled not only to one 
another but to the tank as well. The same is true when we choose a two-photon initial state 
^(0) = \2)(5a(^h- Now we have ui = U2 = l/\/3, = Ui — — 1/\/6, and 

-^e-'""sin(V6At)|l>(a.A + AaJ. (27) 

Here again we are not able to get rid of the entanglement of qubits with the tank. 

5. With this goal in mind wc consider the case when our system "two qubits + tank" 
has an energy equal to the energy of the single photon and can be in a superposition of three 
states {n — 0) : 

*2 = \0)aa(3b, *3 = \0)(3aab, *4 = \l)(3aPb. (28) 

Now it is not enough the energy to excite the first state (8). The total Hamiltonian 
H — Hq + Hint (7) has three eigen-energies (we drop the vacuum energy cut/'^) 

E^^^ = V2X, = -V2X, E^^^ = 0, (29) 

with the corresponding wave functions 

v2 

For the non-symmetric initial state, ^(0) = \0)aaf3b, when there are no photons in the 
tank and the first qubit is in the excited state, the evolution of the wave function of the 
system is described by the expression 

*(() = ii^^^^|^|0>cA - i^i2M^|0),,„„. + -L sin(v/2At)|l>/?A. (31) 

Again, we are not able to obtain the entanglement of the qubits without the intrusion of the 

tank states. 

However, for the symmetric excitation, 

^{0)^\l)Mb, (32) 

when we have a single photon in the tank and both qubits are in the down states, we obtain 
ui — —l/y/2,U2 — l/V^,Us = 0, so that the wave function of the system in the moment t 
will look like 

^{t) = cos{V2Xt)\l) PaPb + ^sin(V2At)|0)^^^^^^. (33) 

v2 



In the moments of time t = {tt/2 + Tmi)/ujji the single photon initially located in the LC- 
circuit produces the pure Bell state of two qubits, 



^a. = ^(-ir|0)^^^^^^^, (34) 

with the completely disentangled tank. Here m = 0, 1, 2, .. and the frequency of Rabi oscilla- 
tions is: 



u;r = V2\ = kJuoTL.n. (35) 



This process resembles the two-slit interference of the single-photon excitation which prop- 
agates into two qubits (slits) simultaneously. 

6. In the case of Mooij type qubit when the persistent current in the qubit loop 
Iq ~ 450 nA, a coupling coefficient k ~ 1/30, a qubit loop inductance Lq ~ 25pH, Lql^ ~ 
6 X 10"^"^ J, and the resonant frequency of LC-circuit J't = ut/^tt = 10^ Hz, or Hut = 
6.3 X 10^^^ J, we have A ~ 0.5 x 10~^^ J, hujR = 1.4 x 10"^^ J, and the Rabi frequency 
Ir = ^^r/^tt = 1.1 X lO^Hz ^ 0.1/2-. Thus, the pure Bell state is formed for the time 
interval 

1 _q 

tBell = — = 2.3 X 10 "sec. 

^JR 

This operation time is much less than the excitation lifetime, Tph, in the LC-circuit with the 
quality factor Q = 1000 : r^^ = 10~^sec, as well as less than the estimated decoherence time 
of the phase qubits . 

In conclusion, in the framework of Jaynes-Cummings model we have investigated the 
protocol for entanglement of phase qubits inductively coupled to a resonant tank. It is shown 
that in the case of simultaneous coupling the single-photon excitation of the tank creates 
the maximally entangled Bell state of qubits with the completely disentangled vacuum state 
of the resonator. The results are applicable to the case of capacitive coupling, as well as to 
the qubit entanglement through optical cavities. 
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